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MODELING  OF  THE  H  +  NF2  REACTIVE  FLOW 

Herbelin  and  Cohen*  proposed  that  an  addition-elimination  rather  than  an 
abstraction  mechanism  might  be  the  preferred  reaction  channel  for  the  produc¬ 
tion  of  electronically  excited  NF(a*A)  from  the  ground-state  reagents  R(»  CH^, 
H  or  D)  and  NF2(2Bj): 

R(2X)  +  NF  (2B  )  ♦  RNF2(2A)  RF  +  NF (a*A  )  (I) 

It  was  found  that  the  electronically  excited  NF(b*E+),  which  is  responsible 
for  the  green  (5288  A)  emission,  can  be  produced  from  the  NF(a*A),  either  by  a 
vibration-electronic  energy-transfer  process  [reaction(2) ]  involving  HF(v  >  2) 
or  an  electronic-electronic  energy-transfer  process  [reaction(3) ]  using 
O^a^  ). 

HF(X1£+,  v  >  2)  +  NFCa1*)  ♦  HF(X1E+,  v  -  2)  +  NF(b1X  +  )  (2) 


02(3^)  +  NF(a*A)  -*■  02(X3E")  +  NF(b1E+)  (3) 

Recently,  Kwok,  Herbelin,  and  Cohen2  measured  the  cross  sections  for  the 
removal  of  NF(b*E+)  by  HF(v  =  0),  H  and  D  atoms,  and  the  radiative  lifetime 
for  the  NF(b*E+)  state.  They  reported  values  for  the  rate  coefficients  for 


*J.  M.  Herbelin  and  N.  Cohen,  "The  Chemical  Production  of  Electronically 
Excited  States  in  the  H/NF2  System,”  Chem.  Phys.  Lett.  20  (6),  (1973), 
p.  605. 

Tl.  A.  Kwok,  J.  M.  Herbelin,  and  N.  Cohen,  "Collisional  Quenching  and 
Radiative  Decay  Studies  of  NF(a*A)  and  NF(b  E+),"  in  Electronic  Transition 
Lasers,  J.  I.  Steinfeld,  ed.,  MIT  Press,  Cambridge,  Massachusetts  (1976), 
pp.  8-12. 
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the  deactivation  of  NF(a*A)  by  the  above  species,  but  were  unable  to  make  a 

direct  determination  of  the  radiative  lifetime  for  this  NF  state.  Since  a 

knowledge  of  this  radiative  cross  section  is  essential  to  the  quantitative 

analysis  of  this  system  and  to  the  determination  of  the  system's  potential  as 

a  chemical  laser,  we  have  attempted  to  extract  this  information  from  the 

q 

previously  reported  NF  density  plots. 

The  number  densities  of  HF  and  NF  were  observed  as  a  function  of  axial 
distance  on  the  large-diameter,  fast-flow,  medium-pressure  flow  tube 
facility. ^  Those  experimental  results  are  summarized  in  Fig.  1.  Since  a 
radiative  decay  rate  of  NF(a^A)  has  not  been  experimentally  established,  we 
report  photon  flux  instead  of  density  alone. 

Our  starting  point  in  analysis  was  a  kinetic  model  developed  in  the 
Aerophysics  Laboratory  that  closely  approximates  the  observed  HFfX1^,  v  =  1, 
2,  3)  vibrational  levels  as  a  function  of  the  flow  distance.^  Also  shown  in 
Fig.  1  is  a  comparison  of  the  theoretical  and  experimental  HF  density  plots^ 
for  a  standard  flow  case.  The  initial  conditions  for  this  computation  are  as 
follows : 


[F] 

=  6. 0  x 

10"11 

mol/cm' 

[h2] 

=  5.9  x 

o 

T 

o 

mol/cm 

[Ar] 

=*  5.3  x 

10~®  mol/cm^ 

^M.  A.  Kwok,  J.  M.  Herbelin,  and  N.  Cohen,  Colllslonal  Quenching  and  Radiative 
Decay  Studies  of  NF(A^A)  and  NF(b^2~t~),  TR-0077(2940)-6,  The  Aerospace 
Corporation,  El  Segundo,  California  (15  April  1977),  p.  16. 

Si.  A.  Kwok  and  R.  L.  Wilkins,  private  communication.  This  is  an  earlier 
version  of  model  discussed  in  R.  L.  Wilkins  and  M.  A.  Kwok,  "Temperature 
Dependence  of  HF(Vi  -  1)  +  HF(v2  "  0)  Vibrational  Relaxation,"  J.  Chem.  Phys. 
70,  1705  (1979). 
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Observed  Number  Densities  as  a  Function  of  Flow  Tube 
Axial  Position  Z 


P  -  1.0  Torr 


T  -  298°K 

U  *  6.94  x  103  cm/sec 

These  conditions  are  determined  experimentally.  The  empirical  velocity  U  is 
1.6  times  the  average  velocity  U  determined  from  the  continuity  equation.  In 
previous  studies,  this  velocity  has  been  found  effective  in  establishing  the 
characteristic  time  "scale"  of  kinetic  phenomena.3  The  kinetic  and  flow  model 
fits  HF(v  «  2)  and  HF(v  =  3)  and  H  densities  well  out  to  30  cm.  The  HF(v  =  1) 
is  consistently  predicted  low,  although  the  relative  decaying  spatial  trend  is 
duplicated.  However,  HF(v  ■  1)  is  not  important  in  the  following  kinetics 
discussions. 

A  nonequilibrium  chemistry  and  constant  pressure  flow  computer  code 
(NEST)3  was  used  to  model  the  full  cross  section  or  radial  elements  of  the 
axisymmetric  cylindrical  tube  flow.  A  second  code  was  used  to  integrate  and 
average  properties  in  the  radial  elements  along  the  tube  diameter  while  cor¬ 
recting  for  the  paraboloidal  velocity  profile.  Corrections  for  radial  varia¬ 
tions  were  found  to  be  small.  Peak  HF(v  *  2,  3)  densities  were  altered  less 
than  15%  of  values  found  by  running  the  NEST  code  for  the  full  cross  section. 
Peak  NF  values  were  reduced  <16%.  In  the  regions  of  decay,  all  densities  were 
adjusted  slightly.  In  Fig.  2,  typical  computed  radial  variations  for  key 
species  are  shown  at  z  =  10  cm  for  the  final  modeled  case.  We  concluded  from 
these  studies  that  species  were  sufficiently  uniform  in  cross  section  in  the 
zone  of  NF(a^A)  and  NF(b^I+)  production.  In  particular,  previous  observa¬ 
tions3  suggested  that,  once  produced,  NF(a^A)  should  be  insensitive  to  wall 
ef  fects . 


5E.  B.  Turner,  G.  Emanual,  and  R.  L.  Wilkins,  The  NEST  Chemistry  Computer 
Program,  TR-0059(6240-20)-l,  The  Aerospace  Corporation,  El  Segundo, 
California  (30  July  1970). 
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In  Fig.  3,  we  attempt  to  fit  the  observed  A(a^A)  [NF(a^A)]  results.  The 

1 3 

rate  coefficient  for  reaction  (1)  with  R  =  H  was  estimated  to  be  2  x  10 
cm^/mol-sec,  based  on  its  similarity  to  the  H  +  NO2  ♦  0H(^JI,  v  =  0)  +  NO 
reaction. ^  Through  the  use  of  this  value  and  the  measured  rate  coefficient 
for  the  removal  of  NF(a^A)  by  H  atoms,  we  were  able  to  reproduce  the  shape  of 
the  experimentally  measured  [NF(a^A)]  curve.  We  chose  an  initial  value  of 
[NF2]  =  1.33  x  10  ^  mol/cmJ  for  the  case,  but  because  the  H  atom  concentra¬ 

tion  was  excessive,  the  shape  is  insensitive  to  this  concentration  up  to 
values  of  10  11  mol/cm  .  As  can  be  readily  seen  from  Fig.  3,  however,  the 
curve  shapes  are  quite  sensitive  to  the  value  of  the  rate  coefficient  for 
reaction  (1).  The  original,  choice  for  [ NF 2 1  yields  the  best  agreement  with 
experiment. 

We  next  introduced  the  transfer  process  (2)  with  a  rate  coefficient  of 
1  ?  3 

k2  =  5  x  10  cm  /mol-sec,  calculated  from  the  measured  value  for  the  reverse 

1  "f*  9 

process,  the  deactivation  of  NFCb^E  )  by  HF(v  =0).  This  implicitly  assumes 
that  electronic-vibrational  energy  transfer  is  the  process  involved.  By  using 
this  value,  we  find  that  the  resulting  computed  zone  plots  have  shapes  similar 
to  the  "A"  curves  in  Fig.  4,  with  a  peak  value  of  [NF(b*£+)]  =  4  x  10-^ 
mol/cm  at  z  =  16.5  cm  downstream.  This  can  be  compared  with  the  experimental 
value  of  8  x  10-1^  mol/cm^  as  shown  in  Fig.  4.  By  varying  the  initial  [NF2] , 
it  was  possible  to  vary  the  computed  peak  [NF(a^A)]  and  [NF(b*£+)]  values  over 
an  order  of  magnitude  range  without  appreciably  affecting  the  shapes  of  the 
curves  for  NF,  HF,  or  H. 

The  shape  of  the  experimental  curve  for  NF(b*E+)  relative  to  the  computed 
curve  "A”  indicated  that  the  following  additional  transfer  process  between  HF 
and  NF  is  also  important: 

HF(X1E+,  v  =  3)  +  NF(a^ A)  *  HF(X1£+,  v  “  1)  =  NF(bJE+)(4) 


^D.  L.  Baulch,  D.  D.  Drysdale,  D.  G.  Horne,  and  A.  C.  Lloyd,  Evaluated  Kinetic 
Data  in  High  Temperature  Reactions,  2,  London;  Butterworths  (1973)  p.  445. 


Model  Fit  of  H  +  NF2  Production  of  NFCa1^)  Showing 
Sensitivity  of  Rate  Constant 


NUMBER  DENSITY,  (mol-cm 


Previous  experience  with  other  transfer  processes  between  NF(a*A)  and  HF, 
02(a^A)  and  more  recently  I*  atoms2  indicates  that  energy  resonance  could  be 
used  to  estimate  a  value  for  the  rate  coefficient  for  reaction  (4).  The 
exothermicity  for  process  (2)  is  298  cm”^  as  written,  assuming  all  species  to 
be  in  the  J  =  0  rotational  state.  Transfer  process  (4)  is  50  cm-1  exothermic. 
Since  this  is  less  than  the  110  cm  ^  for  the  transfer  process  between  1  atoms 
and  NF,  for  which  a  gas  kinetic  cross  section  has  been  measured,  we  assumed 
that  the  vibrational  energy  of  the  HF(v  =  3)  would  transfer  in  one  out  of 
every  two  collisions.  By  using  an  approximate  collisional  diameter  of  3  A  for 
all  species  involved,  we  arrived  at  a  room  temperature  rate  coefficient  of 

1  O  Q 

k(4)  =  6  x  10IJ  cmJ/mol-sec  in  the  exothermic  direction.  The  endothermic 
reverse  reaction  rate  coefficient  was  then  calculated  to  be  k^_^  =  6  x  10^ 
exp(-142/RT). 

Inclusion  of  this  transfer  process  into  the  model  resulted  in  considerably 
improved  agreement  between  the  shapes  of  the  computed  and  measured  curves.  By 
trial  and  error,  we  adjusted  this  rate  coefficient  and  the  initial  NF 2  concen¬ 
tration  to  the  final  values  of  =  2.6  x  10*2T122,  k_^  =  2.6  x  lO^T1^2, 
exp(-142/RT)  cm^/mol-sec,  and  [NF2]  =  4.2  x  10~^  mol/cm^  to  arrive  at  the 
agreement  between  the  computed  curves  B  and  the  experiment.  Only  a  50%  adjust¬ 
ment  from  the  initial  estimate  was  required  to  achieve  this  excellent  fit. 

Calculations  were  also  performed  to  explore  the  effect  of  the  additional 
transfer  process, 

HF(X!Z+,  v  =  2)  +  NF(b1Z+)  t  HF(X1I+,  v  =  4)  +  NF(a*A  )  (5) 
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J.  M.  Herbelin  and  M.  A.  Kwok,  Electronic  Angular  Momentum  Transfer  in 

+  NF(a^A)  System,  TR-0076(6240-40)-4,  The  Aerospace  Corporation, 
El  Segundo,  California  (17  March  1976). 
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We  found  that  the  computed  NF(b^E  +  )  profiles  were  quite  insensitive  to  the 
process  even  when  a  gas  kinetic  rate  coefficient  was  assumed.  Therefore,  no 
conclusion  concerning  this  process  is  possible. 

When  this  agreement  between  observed  and  computed  values  of  [NF(b*£+)] 

and  its  axial  variation  has  been  obtained,  the  computed  population  of  the 

N'F(a^A)  can  be  used  to  infer  the  radiative  lifetime  (inverse  of  the  radiative 

decay  rate)  of  the  NF(a^A)  NF(X^E  )  transition.  We  simply  take  the  ratio  of 

9  1  +  1 

the  observed  integrated  band  intensities  from  the  (b  £  )  and  (a  A)  states 
which  is  proportional  to  the  ratio  of  wavelengths,  populations,  and  radiative 
decay  rates. 


l(  b1E+  ->•  X3I  )  exptl  _  \(  a*A]  Nf  b1!*)  A[  b1!*) 
l(  a*A  X3Z~)  exptl  x(b1E+)  N(alA)  a(  a  'a  ) 


(6) 


We  find  that  the  ratio  of  radiative  decay  rates 


A-^aA-|  =  0.015 

A(  blZ+) 

which  yields  a  value  of  A(a^A)  =  0.94  sec  ^  when  the  value  of  A(b*£+)  =  67 
-1  2 

sec  as  measured  by  Kwok,  et  al.  is  used.  This  method  requires  only  a 
relative  intensity  calibration.  The  value  for  A(a^A)  was  also  computed 
independent lv  on  the  basis  of  the  absolute  (a-X)  emission  band  intensity. 
Absolute  intensitives  yield  results  such  as  the  photon  flux  plot  in  Fig.  1. 
if  NF(a^A)  is  known,  A(a*A)  is  readily  determined.  In  this  case  A(a*A)  is 
deduced  to  be  1.75  see'*.  The  two  calculated  values  agree  to  within  a  factor 
of  two. 

As  a  final  check  of  results  from  this  analytical  procedure,  we  performed 
an  experiment  in  which  we  went  from  the  hydrogen-rich  regime,  [H]  »  [NF2I,  to 
one  in  which  the  NF2  was  in  slight  excess.  We  accomplished  this  by  decreasing 
the  NF3  flow  and  discharge  power  until  we  reach  a  fluorine  atom  flow,  measured 
by  our  standard  testing,  of  approximately  0.05  of  the  test  case.  The  total 
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pressure  was  maintained  at  the  1  Torr  value  and  excess  deuterium  was  used  to 
produce  a  concentration  of  D  atoms  of  2.3  x  10  mol/cm  .  We  then  introduced 
NF2,  produced  by  passing  N2F4  through  an  oven  prior  to  injection  into  the  flow 
tube.  The  NF2  was  injected  approximately  20  tube  diameters  (2  m)  upstream  of 
the  N2  injector  to  ensure  adequate  mixing.  The  NF2  flow  was  calibrated  by 
supplying  an  excess  of  N  atoms  from  an  additional  microwave  discharge  through 
N2  and  then  measuring  by  the  standard  method  the  F  atoms  produced  by  the 
reaction 


N  +  NF2  +  +  2F  (7) 

The  reduced  concentration  of  the  D  atoms  resulted  in  greatly  extended 
NF(a*A)  zones,  in  which  the  concentrations  of  NF(a^A)  were  approximately  flat 
at  distances  in  excess  of  40  cm.  Figure  5  shows  the  NF(a^A)  density  plot  for 
such  a  NF2  flow  condition.  The  peak  is  now  reached  much  later  in  Z  because 
the  production  rate  kj[NF2)  is  much  slower  than  the  previous  kj[H]  of  the 
[H]  »  [NF.j  regime. 

In  Fig.  6,  we  show  the  variation  of  [NF(a*A)]  deduced  from  the  measured 
peak  signals  with  the  use  of  the  above  deduced  1.75  sec  *  value  for  the 
radiative  decay  rate  of  NF(a^A).  Excess  O2  was  added  in  one  run  to  remove  any 
interfering  N2  first  positive  emission,  which  is  seen  to  have  a  20%  effect. 
Because  of  the  limits  of  our  N2F4  heating  facility,  we  were  not  able  to  pro¬ 
duce  more  than  a  7.5  umol/sec  flow  of  NF2.  We  see  that  the  [NF(a^A>]  rises 
and  appears  to  reach  a  limit  of  3.0  x  10  raol-cm  .  This  value  is  in  the 

mid-range  of  limits  of  H  atom  density.  The  lower  H  limit  was  determined  by  an 
F-atom  titration  and  by  the  assumption  of  full  conversion  to  H  atoms  by  the  F 
+  H2  reaction.  The  upper  H  limit  was  given  by  the  assumed  total  conversion  to 
F  by  the  NF3  passed  through  the  discharge.  It  appears  that  reaction  (1)  is 
the  key  mechanism  in  producing  NF(a^A)  since  its  maximum  density  equals  that 
of  the  limiting  reagent  H  within  a  factor  of  two. 
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[H],  UPPER  LIMIT 


[NFIa^AI],  Z=65cm 
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Peak  NF(a^A)  Density  as  a  Function  of  NF2  Flow  Rate 


This  result  supports  the  above  analyses  as  well  as  the  contention  that 
•eaction  (1)  is  close  to  100%  efficient  in  producing  electronically  excited 
IF(a^A).  Further  work  is  required  to  confirm  this  latter  conclusion. 


APPENDIX 


Since  this  analytical  work  was  done,  independent  experimental  work  has 

8  Q 

been  reported  by  Cheah,  Clyne,  and  Whitefield,  and  by  Malins  and  Setzer. 
The  experimentally  determined  rate  coefficients  (at  room  temperature)  for 
reaction  (1)  can  be  compared  with  our  analytical  deduction  in  the  following 
tabulation: 


10"13  kj (cm3-mol-^sec~^ ) 

This  work  2 

Cheah,  Clyne,  Whitefield  0.90 

Malins,  Setzer  >0.37 

We  consider  our  deduction  to  be  in  good  agreement  with  the  measurement  of 
Cheah,  et  al.,  which  relied  on  direct  observation  by  UV  resonance  absorption 
of  fluorescence  to  monitor  the  disappearance  of  H  atoms  caused  by  the  presence 
of  NF2  radicals.  They  did  not  resolve  the  problem  of  whether  product  state 
NP(a^A)  is  in  a  totally  inverted  condition.  If  this  assumption  is  made,  then 
Cheah,  et  al.  have  succeeded  in  determining  the  reaction  rate  coefficient  of 
reaction  (1). 

The  measurements  of  Malins  and  Setzer  confirm  the  exciting  prospect  that 
over  90%  of  the  H  +  NF2  reaction  produces  the  excited  electronic  state, 
NF(a3A).  However,  their  rate  measurement  depends  on  a  method  of  comparison 


®C.  T.  Cheah,  M.  A.  A.  Clyne,  and  P.  D.  Whitefield,  "Reactions  Forming 
Electronically-excited  Free  Radicals,"  J.  Chem.  Soc.  Faraday  II  76  (1980) 
pp.  711-728. 

*R.  J.  Malins  and  D.  W.  Setzer,  "Rate  Constant  and  Branching  Ratio  for 
Formation  of  Electronically  Excited  NF  by  Reaction  of  H  Atoms  with  NF  by 
Reaction  of  H  Atoms  with  NF2«"  (To  be  published).  We  thank  Dr.  Malins  for 
advance  preprints  of  his  work  and  for  many  helpful  discussions. 
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with  another  "baseline"  reaction,  and  it  may  not  be  surprising  that  their 
absolute  numbers  are  in  only  fair  agreement  with  Cheah,  et  al.  and  this  work. 
Their  method  of  diagnosis  prevents  them  from  observing  all  the  product  states 
of  HF(v),  namely  HF(v  =  0),  which  is  expected  to  be  a  significant  product. 
Therefore,  they  indicate  their  rate  coefficient  could  be  larger. 

Malins  and  Setzer  have  also  attempted  the  first  measurement  of  the 
NF(a^A)  radiative  lifetime  with  a  result  of  5  to  6  sec.  Our  deduced  results 
yielded  0.6  to  1  sec  by  two  different  approaches.  We  estimate  the  Malins 
number  to  be  more  reliable,  but  the  factor  of  5  disagreement  suggests  that 
even  more  experimental  work  should  be  done  on  this  interesting  reaction 
system. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  Is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laborato¬ 
ry  personnel  in  dealing  with  the  many  problems  encountered  in  the  Nation's 
rapidly  developing  space  systems.  Expertise  in  the  latest  scientific  develop¬ 
ments  is  vital  to  the  accompl 1 shment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are:  f 

Aerophystcs  Laboratory:  Aerodynamics;  fluid  dynamics;  plasmadynami cs; 
chemical  kinetics;  engineering  mechanics;  flight  dynamics;  heat  transfer; 
high-power  gas  lasers,  continuous  and  pulsed,  1R,  visible,  UV;  laser  physics; 
laser  resonator  optics;  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  re." 't Ions  and  optical  back- 
grounds;  radiative  transfer  and  atmospheric  transmission;  thermal  and  state- 
specific  reaction  rates  in  rocket  plumes;  chemical  thermodynamics  and  propul¬ 
sion  chemistry;  laser  isotope  separation;  chemistry  and  physics  of  particles; 
space  environmental  and  contamination  effects  on  spacecraft  materials;  lubrica¬ 
tion;  surface  chemistry  of  Insulators  and  conductors;  cathode  materials;  sen¬ 
sor  materials  and  sensor  optics;  applied  laser  spectroscopy;  atomic  frequency 
standards |  pollution  and  toxic  materials  monitoring. 

Electronics  Research  Laboratory:  Electromagnetic  theory  and  propagation 
phenomena;  microwave  and  semiconductor  devices  and  Integrated  circuits;  quan¬ 
tum  electronics,  lasers,  and  electro-optics;  communication  sciences,  applied 
electronics,  superconducting  and  electronic  device  physics i  mlllloeter-wave 
and  far- infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  composite 
materials;  graphite  and  ceramics;  polymeric  materials;  weapons  effects  and 
hardened  materials;  materials  for  electronic  devices;  dimensionally  stable 
materials;  chemical  and  structural  analyses;  stress  corrosion;  fatigue  of 
metals. 

Space  Sciences  Laboratory:  Atmospheric  end  Ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  alrglov;  magnetoapherlc  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  In  the  magnetosphere;  solar  physics,  x-ray  astronomy;  the  effects 
of  nuclear  explosions,  magnetic  storms,  and  solar  activity  on  the  earth's 
atmosphere,  Ionosphere,  and  magnetosphere;  the  effects  of  optical,  electromag¬ 
netic,  and  particulate  radiations  In  space  on  space  systems. 


